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Abstract 

A next-to-leading order QCD analysis of deep inelastic scattering data is performed al- 
lowing for contributions due to a light gluino. We obtain the values of as{M'^)±6af"'^ = 
0.108±0.002, 0.124±0.001, 0.145±0.009 for QCD, SUSY QCD with a Majorana gluino 
and a Dirac gluino respectively. The value of as{M'^) obtained in SUSY QCD with a 
Majorana gluino best agrees with the direct measurements of as{M'^) at LEP. 



The most precise measurements of asiQ"^) are provided by deep inelastic lepton-hadron 
scattering and e^e~ experiments and extend over many generations of scale providing a 
test of QCD [^ . The results of the QCD analyses performed in the different deep inelastic 
scattering experiments may be expressed in terms of a value for as(M|) (by a next-to- 
leading order (NLO) relation, for example) and thus directly compared with the results 
obtained in e+e~ experiments at LEP. Different values of UsiM"^) are obtained from deep 
inelastic and e+e~ experiments. The question has arisen whether this difference might be 
due to the presence of a light gluino B R] for which windows remain open in the range 



m 



a 



< 



5 GeV i IQ. 



Unlike direct measurements of asiM"^) in LEP experiments, the QCD analysis of deep 
inelastic scattering data requires an assumption on the /3-function to solve the evolution 
equations. As shown in [^] the re-evaluation of as(M|) resulting from the DIS-data under 
the assumption of the NLO /3-function in the presence of a light Majorana gluino yields 
a better agreement with data. 

However, the effect of a light gluino in the mass range m- ~ 5 GeV on the scaling 
violations of structure functions is not only due to a modification of the /5-function, but also 
due to modified evolution equations containing even different parton densities. Since the 
allowed gluino mass is rather small, a thorough study of this effect requires a complete NLO 
QCD analysis of the deep inelastic data including light gluinos|. This is the purpose of the 
present paper. 



The Running Coupling Constant 

Before we study the evolution of deep inelastic structure functions we summarize the effect 
of a gluino on as{Q^)- The renormalization group equation for a^ 
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yields the solution 
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with the function^ ^^"^'(x; Pi) having the form in next-to-next to leading order (NNLO) 
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^Different possibilities for a dedicated search for light gluinos have been also proposed in [0, 0, R| recently. 
Light gluinos are also favored in some approaches in string theory (see [pj). 

^The effect of a gluino with m~ ~ 5 GeV on data expected in the high Q^ range at HERA, decoupling at 
low Q"^ range, was studied in in LO recently. 

■^Because the explicit expression for (|^) depends on the size of /3o,i.2.... ^^"H^'i f^i) i^ understood to be a 
complex valued function. 



In the limit /52 — > (NLO) one obtains 
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where we choose C = following [03. The coefficients /3i are given by [|2 
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(A^'j, A'^) describes the gluino contribution to /?2- ^/ denotes the number of active quark 
flavours and N~ = 0, 1,2 refers to QCD, supersymmetric QCD with Majorana gluinos, and 
Dirac gluinos, respectively. Eq. (0) can be solved iteratively. If as(Qo) i^ chosen by q;s(M|) 
one can calculate as{Q'^) as a function of /3o and /3i, and the dependence on A is implicit 
only through ^^(Qg). The A dependence can be made explicit through the definition 
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in the MS scheme in NLO0. Eq. (^ relates the values of A and a^ at a given scale directly. 

When quark or gluino mass-thresholds /i ~ 2mj are passed, as(Q^) is kept continuous 
as the values of /3j change. To be more precise, running mass effects should be accounted 
for when passing quark or gluino mass thresholds |T^. We will neglect these effects in the 
present analysis as smallQ 

In figure 1 a comparison of the NLO solutions @ using q;s(Qo) — <^s{Mz) = 0.122 for 
A'- = 0, 1, 2, m~ = 0, 3, 5 GeV, and different measurements of a^ is shown. Here we assumed 
rric = 1.6 GeV and nib = 4.75 GeV for the values of the heavy quark masses. For the case 
of QCD also the curve obtained in NNLO is given. That a light Majorana gluino is an 
apparently better description of as{Q^ > 16 GeV"^) than QCD has been previously observed 
in 1^. Calculating confidence levels for the depicted experimental data in the complete Q"^ 
rangeQfrom the x^ values for the three different hypothesis Ng = 0, 1, 2, however, one obtains 
78 % CL for QCD, 50-83 % CL for the case of a Majorana gluinos (0 < m^ < 5 GeV), and 
< 0.5 % CL for the case of a Dirac gluino. This estimate, however, assumes, that the a^ 
values have been extracted in a way which is essentially insensitive to the value of Ng, since 
most of these values were obtained using QCD matrix elements as an input. 



QCD Analysis 

Scaling violations of the structure functions measured in different deep inelastic scattering 

■* This choice is often caUed definition of as{Q^) in the MS scheme, ahhough /3o,i are scheme independent 
and only /32 eq. M) refers to a MS resuh. Other possible choices are discussed in pi[ . 

^A corresponding relation in NNLO follows from (0,0) directly. 

^ As has been shown recently 111] this modification leads to a shift Aas(M|) ^ 0.001 unless a quark 
mass threshold coincides with that of the gluino. 

^The as value determined from cc IP — IS" mass splitting using a lattice calculation |15[| was excluded. 



experiments can be described^ by the following set of NLO evolution equations including the 
effect of a light gluino, 
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The gluino density g{t, x) emerges in the singlet evolution equation. In the description of the 
singlet and non-singlet distributions we assume that the squark contributions decouple p 
due to the current experimental limit rriq > 74 GeV [Q . The functions 
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are expressed by the splitting functions P'^^\x) up to NLO. Note, that the leading order 
matrix elements P ~ ix) and P~ ix) of P^^'ix) vanish. 

qg ^ ' gq ^ ' ^ ' 

The different combinations of quark densities in (|^) are 

(it\t,x) = qi{t,x)±qi{t,x) 
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The evolution variable is t = — (2//3o) loglasiQ"^) / as{Ql)], with Qq defining the scale at which 
the input distributions g| (x, Ql), G{x, Qq), and g{x, Ql) are parametrized. describes the 
Mellin convolution. 



The splitting functions P/ (x) in (0,|^) were given in |T8|, |T9| where in NLO the DR scheme 
was used. In the present analysis these results were translated into the MS scheme. 

The QCD analysis was performed with a modified version of the CTEQ parton distri- 
bution evolution program [^]. The input shapes for the parton distributions at Q^ were 
chosen to be 
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^The data analysis performed on this basis will include the data currently measured at HERA ||l^ also. 
Although these are small x data they also exhibit logarithmic scaling violations. 



x{d — u){x) = —x'^~{l — x)^'{l + 5-\^ + 'j-x) 

xs{x) = Asx"" {1 - x)'^' 

xG{x) = Agx°«(1-x)^^(1 + 7gx) 

xg{x) = AgX^^il-xfs (ir 



The sensitivity of the parameters in { pA\j to the fit resuhs varies, and some of them may 
be fixed or can be related to one another. The 10-parameter parametrization used in our 
analysis was the result of a study to minimize the number of shape parameters needed to 
fit the current set of experimental inclusive data in a global QCD analysis. This is desirable 
as fewer shape parameters are clearer to analyze and reduce spurious correlations between 
themselves. We set a„ = a^ = 0.5, j3u = 4.0, j3d = 3.0, /3+ = /3_ = /5g + 1, a+ = ac, and 
7_,_ = 'fQ. The sum rules J^ u^{x)dx = 2, ^^ xdy{x)dx = 1, and momentum conservation 
reduce by three the number of input fitting parameters. The parameters of the strange 



quark distribution have been determined from the CCFR dimuon data |21[] as Ag = 0.114, 
as = —0.114, and Ps = 6.87 and are treated as fixed. Only for m- < Qo = 1.6 GeV 
do we have an initial three parameter nonzero gluino distribution. Otherwise the gluino 
distribution was radiatively generated as were the heavy flavor distributions. 

We allowed the relative normalizations of the different experimental data sets (seven 
parameters) to float under the fit subject to suitable constraints reflecting experimental 
normalization uncertainty in the mean. 

A and the parameters of eq. (|Tl|) were determined by a MINUIT [^ minimization of a x^^ 
distribution using the following data sets: SLAG F^^, Ff, Ff/F^^, NMC F^, F^^, F^'^/F^'^ 
for E^, = 90 and 280 GeV, BGDMS F^'^, F^^, GGFR F^, xF^, and ZEUS and HI F^p [g|, 
as input. The statistical and systematical experimental errors were added in quadrature, 
and relative systematic effects between the different experiments were accounted for by 
normalization factors determined in the fit. For a given gluino mass hypothesis it was 
demanded, that W^ = Q'^{1 — x)/x > 4m? leading to a constraint on the data used. For 
the comparison of different possible hypotheses we considered the cases N~ = (QGD), and 
A'- = 1,2 with m~ = 0, 3, and 5 GeV. Table 1 summarizes the results on A and a;s(M|) and 
table 2 contains the shape parameters determined for some characteristic cases and the fitted 
relative normalizations of the different data sets used in the analysis. The latter parameters 
values differ by only up to 5 % from unity. Note, that the values of Ai^^ (the value above 
the 6-quark threshold) and q;s(M|) given in table 1 are directly related by eq. (H) and no 
reference to a "typical" value of (Q^) characterizing the analysed data sets is needed]^. 



Discussion 

The values of x^/NDF obtained in the analysis of all the deep inelastic scattering data do not 
differ significantly (see table 1). Thus, the analysis does not yield an a priori preference to one 
of the cases N~ = 0,1,2. The fitted values of A^ are expressed in terms of as{M'^) in table 1. 
as(M|) allows a more direct comparison of the different results, because Mz ^ rrii^mg. 
Besides the quoted statistical error (5af"*, there is a theory error due to the assumption on 



^Such an assumption would be required for an illustration of the obtained result in the way shown in 
figure 1, but will not be given due to the uncertainty in the abscissa. 



the factorization scale [^], (5as(M|)*'*^°^, which has been estimated to be ~ 0.005. This error 
dominates the statistical error determined by the analysis. From direct measurements at LEP 
one obtains as(M|) = 0.122 ±0.006 [lb]. Comparing the results on as(M|) given in table 1 
with the LEP result, we derive the confidence levels of 8 %, 70 %, and 5 % for the hypotheses 
Ng = 0, 1, and 2, respectively, for the agreement of both values. Here 6af''^°^ = 0.005 was 
added to the corresponding statistical errors in quadrature in the definition of x^- 

The resulting up quark, gluon and gluino distributions at Q^ = 10 GeV^ and 100 GeV^ 
of QCD, SUSY QCD with a massless, 3 and 5 GeV gluino are shown in Fig. 2a-f. The 
distributions ususy{x) and gsusrix) are softer at Q^ < 10 GeV^, the faster running coupling 
of QCD compensating for the initial depletion of larger x distributions in the SUSY case. 
At small X, X < .01, the sea and gluon distributions in the SUSY case are enhanced by the 
splitting of larger x gluinos. 

The fit results are not very sensitive to the parameters of the gluino distribution and allow 
large distributions at small x where the data is not yet very precise however. The gluino 
distribution is forced by the evolution equations to behave like a sea quark, the existing 
DIS data forces it to be quite soft. The magnitude of the gluino distribution is strongly 
dependent on rrig over a large range of Q^, the dependence decreases logarithmically with 
increasing scale. 

If Majorana gluinos exist in the mass range rrig < 5 GeV, they can be searched for in single 
and pair production in deep inelastic ep-scattering and photoproduction at HERA. The 
process of single gluino production will moreover allow the derivation of direct constraints 
on the gluino distribution. 
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Table 1: Results of the QCD analysis for A and as(M|) 
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relative normalizations 


SLAG 

NMC v^ = 90 GeV 

NMC Vs = 280 GeV 
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GGFR 

ZEUS 

HI 


1.002 
1.004 
1.016 
0.984 
0.971 
1.010 
0.961 


0.998 
1.003 
1.010 
0.978 
0.966 
1.002 
0.952 


0.999 
1.003 
1.013 
0.981 
0.967 
1.004 
0.956 



Table 2: Results of the QGD analysis for the parameters of the input distributions at 
Ql = 2.56 GeV=^ 
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Figure 1: Comparison of different theoretical predictions for UgiQ'^) with experimental results of 
as [||]. The full curves denote the NLO solution of eq. (§) for N~ = 0,1,2 with m~ = taking 
as{Ql) = OLs{Ml) = 0.122. The dash-dotted hue denotes the NNLO solution in the case of QCD. 
The dashed and dotted lines describe the cases m~ = 3 and 5 GeV, respectively. 
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Figure 2: Fit results for the parton distributions, a) xu{x,Q'^), b) xG{x,Q'^), c) g{x,Q'^) for 
Q2 = lOGeV^. The thick sohd, thin sohd and dashed hues denote QCD, SUSY QCD with a 
massless Majorana gluino and SUSY QCD with a Majorana gluino of nig = 3 GeV respectively. 
d)-f) are same as a)-c) with Q^ = 100 GeV^. The dotted line denotes SUSY QCD with a Majorana 
gluino of nig = 5 GeV. 
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